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Abstract 

Inspired by recent lattice measurements, we determine the short-distance (a <C r <C 1/itT) 
as well as large- frequency (1/a S> 00 ^> irT) asymptotics of scalar (trace anomaly) and 
pseudoscalar (topological charge density) correlators at 2-loop order in hot Yang-Mills theory. 
The results are expressed in the form of an Operator Product Expansion. We confirm and 
refine the determination of a number of Wilson coefficients; however some discrepancies with 
recent literature are detected as well, and employing the correct values might help, on the 
qualitative level, to understand some of the features observed in the lattice measurements. 
On the other hand, the Wilson coefficients show slow convergence and it appears uncertain 
whether this approach can lead to quantitative comparisons with lattice data. Nevertheless, 
as we outline, our general results might serve as theoretical starting points for a number of 
perhaps phenomenologically more successful lines of investigation. 
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1. Introduction 



Many of the most interesting physical properties of a finite-temperature system are of an 
"infrared" type, i.e. to be extracted from the long-distance or short-frequency limit of ap- 
propriate 2-point correlation functions. For instance, in a weakly coupled Yang-Mills theory 
with the gauge coupling g, spatial correlation lengths originate at length scales r > l/(gT) or 
r>Tr/(g 2 T), depending on the global quantum numbers of the operator considered pQ. At 
the same time, transport coefficients, which reflect the real-time response of the system to 
small perturbations, arise at frequencies w<g 4 T/7r 3 (cf., e.g., ref. |2|). Infrared observables 
may either be genuinely non-perturbative [3j|4], or they do possess a weak-coupling series up 
to some order, but it is slowly convergent at temperatures relevant for heavy ion collision ex- 
periments (cf., e.g., ref. [5]). Hence, perhaps with a few exceptions (cf., e.g., ref. [6j), infrared 
observables need eventually to be determined via non-perturbative lattice simulations. 

Despite the stated general picture, circumstances exist as well under which "ultraviolet" 
observables, measured at short distances (r <1/ttT) or large frequencies (uj>ttT), are of 
physical interest. As an example on the former case, we may mention that in heavy quarko- 
nium physics, the system has an additional "external" scale, the heavy quark mass, M. Given 
that normally ttT <C M, the inverse Bohr radius 1 ~ a s M <C M of a quarkonium state 
could well be of the same order as the temperature, r^ 1 ~ irT. Changes in quarkonium 
properties caused by a finite temperature could therefore be due to thermal modifications of 
the quark-antiquark potential at r ~ r B ~ l/(irT), in which case the weak-coupling expan- 
sion may converge faster than at large distances r > l/(gT) [7\. As an example on the latter 
case, we remark that all lattice estimates of a spectral function, p(oj), from whose intercept, 
lirmj_>o p(oj)/oj, transport coefficients are determined, rely on "inverting" the relation 



where (3 = 1/T and G(f), < f < 1, is a measured Euclidean correlator along the compact 
time direction. Since the "kernel" multiplying p(oj) in eq. (jl.ip only depends on oj through 
/3uj, it is clear that p{oj) in the regime oj ~ irT gives an important contribution to G(f), and 
needs to be well understood before infrared sensitive contributions from the range oj <C ttT 
can be reliably extracted. 

In the present note we study the ultraviolet regime of certain 2-point correlation func- 
tions, and even take it to its extreme limit: not only do we consider r <l/(irT) but in fact 
r <C 1/(ttT); not only oj > ttT but in fact oj S> vrT. Only the ultraviolet cutoff, e.g. the lattice 
spacing a, is assumed to be even farther in the ultraviolet than our physical scales. In this 
situation the correlation functions can be determined within a framework similar to the Op- 
erator Product Expansion [8], as has recently been discussed in the Euclidean domain [9] and 
also more generally [10] . Apart from theoretical considerations, we would in principle also like 
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to make contact with the Euclidean lattice simulations described in ref. [II]. It is our expe- 
rience, both from Euclidean [7J and, after analytic continuation, Minkowskian [12] domains 
that, when pursued to a sufficient order, the (continuum) weak-coupling expansion could 
work relatively well in the ultraviolet regime; however, the issue needs to be re-investigated 
when a new correlation function is considered or when the computation is organized as an 
Operator Product Expansion, and these are some of the goals of the present study. 

The plan of this note is the following. In sec. [2] we define the basic observables considered. 
Section [3] contains an outline of the method used; the general results are given in sec.|U The 
case of short distances is discussed more specifically in sec. El and that of large frequencies 
in sec. EJ We also briefly comment on the relation of our work to recently discussed "sum 
rules" in sec. [7j Section [8] offers a summary and outlook, whereas the three appendices collect 
together some details needed in the main text: the definitions and asymptotic expansions of 
the "master" sum-integrals appearing in pure Yang-Mills theory are listed in appendix A; a 
number of thermodynamic potentials playing a role in our study are given in appendix B; 
and fermionic effects are briefly discussed in appendix C. 



2. Setup 

Employing the convention D p = d p — ig B Ap[T a , with T a hermitean generators of SU(iV c ) 
normalized as Tr [T a T b ] = \5 ab , and defining F* v = (2i/c/ B )Tr {T a [D^, D u ]} = d^A* - 
d v Ap\ -\- g B j ahc A h pA c v , the dimensionally regularized Euclidean action relevant for pure Yang- 
Mills theory at a finite temperature T = 1/(3 reads 



S E = I dr ld 3 - 2 ^\ 7 F^ u \ . (2.1) 



The trace of the corresponding energy-momentum tensor is ^£-^F®F®, where D = 4 — 
2e is the space-time dimensionality. In order to define the operators whose correlation 
functions we are interested in, we note that it is the "geometric" structure g\F^ v Fp a = 
— 2Tr {[Dp, D u ][D p , D a ]} which requires no renormalization at the order of our computation. 
So, we define the gauge invariant scalar and pseudoscalar operators 

e = c e g 2 B F^F^, X = c x e^ pa glF^F^, (2.2) 

where the D-dimensional Euclidean e^ vp(T is handled as in ref. |13j (as long as the procedure 
is self-consistent the precise regularization has no impact on the final results). In analogy 
with ref. [11] the coefficients are chosen as 

£-4 b hg 2 
° = i^ = -7- — + -' (2 - 3) 

C ^ 64^' M 
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but in most of what follows we do not need to specify their values. For eq. (|2.3p . we have 
written 



2 2 2e 

9b = 9 r/i 



bog 2 , (b\ h \ 4 



1 -^- + Uf-si' + - 



(2.5) 



where 



and g 2 is the dimensionless renormalized coupling constant, evaluated at the MS scheme 
renormalization scale p, (fi 2 = /2 2 e 7E /47r). 

With this notation, the Euclidean correlators considered are defined as 

G e {x) = (0{x)0(0)) c , G x (x) = (x(^)x(O)) , (2.7) 

where (...) c denotes the connected part, and the expectation value is taken at a finite tem- 
perature T. (The disconnected part of Gq is (0) 2 , where (6) = e — 3p is trace of the energy- 
momentum tensor; in a general regularization scheme (9) is ultraviolet divergent. In the 
following we need the finite thermal part thereof, denoted by (e — 3p)(T).) We also consider 
the corresponding Fourier transforms, 

G e {P) = [ e- W x G e {x) , G X (P) = [ e~ ip - x G x {x) , (2.8) 

J X J X 

where short-distance singularities are regulated dimensionallyQ Finally we denote 

AG e (P) = G 6 {P) - Gj=°(P) , AG X (P) = G X (P) - G T X =\P) , (2.9) 

subtracting the zero-temperature parts at a fixed P. In Euclidean signature, the components 
of the four-momentum and spacetime coordinate are denoted by 

P = (Pn,p), P=|p|, x = (t,x), r = |x| , (2.10) 
where p n = 2i:Tn are bosonic Matsubara frequencies. 



3. Method 

Carrying out the Wick contractions in the 1-loop and 2-loop Feynman graphs contributing to 
eq. (I2.7P (cf. ref. [13] or fig. [T]) and using the notation of appendix A for the various "master" 

^Tor a careful general analysis of short-distance singularities in the case of G x , see ref. [14] . 
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Figure 1: The graphs contributing to the connected 2-point correlators defined in eq. (|2.7p . up to 
2-loop order. The wiggly lines denote gluons; the small dots the operators 8 or x (cf- eq. (|2.2p ); and 
the grey blob the 1-loop gauge field self-energy. Graphs obtained with trivial "reflections" from those 
shown have been omitted from the figure. 



sum-integrals appearing in the result, we find the bare expressions (d A = N 2 
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(3.1) 



(3.2) 



Here, for brevity, structures containing $ Q 1, which vanishes exactly in dimensional regular- 
ization, have been omitted (£ Q denotes a sum-integral with bosonic Matsubara frequencies). 

The goal now is to obtain asymptotic expansions for the functions in eqs. (|3.ip . (|3.2p . valid 
in the Euclidean domain P 2 = Pn+p 2 (ttT) 2 . (Interestingly, as has recently been reviewed 
in ref. |15| . similar expansions may play a role also in the computation of high-order vacuum 
graphs at finite temperature.) 

In order to obtain the expansions, we first carry out all the Matsubara sums, which can 
be done exactly. This has the effect of setting one or two of the propagators "on-shell"; 
the corresponding line is weighed by the Bose distribution. The coefficient can be identified 
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as a zero-temperature amplitude or 1-loop integral, with special kinematicsl§ For instance, 
introducing the notation 

[•••]q = o ^2 {•••}> 1"]q,r=J ^2 {•••}) 



(3.3) 

q n =±iq q n =±iqr n =±ir 

it is straightforward to verify that the sum- integral 1^ (cf. eq. (|A.10p ) can be re-expressed as 

p4 



Ic 
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(3.4) 



Q,R 



_(Q-R) 2 (R-P) 2 (Q-R) 2 (Q-R-P) 2 (R - P) 2 (Q + R - P) 2 _ 

where Jq = J d D Q/(2ir) D and f = J d D_1 q/(27r) £)_1 . It is important to realize that within 
the square brackets we can set Q 2 = for any D, and that therefore anything proportional 
to D-dependent powers of Q 2 vanishes exactly (this is the case particularly when the R- 
integration factorizes from the P-dependence). 

In order to handle the remaining structures, we note that an integration variable appearing 
inside the Bose distribution is always ultraviolet safe (cut off by the temperature). Therefore, 
we can expand propagators as 



1 




[(Q~R) 2 _ 


Q 



1 2Q-R 4(Q • R) 2 
R 2+ ~R 4 ~ + ^ 



(3.5) 



where we also made use of the on-shell condition Q 2 = 0. Such an expansion leads to the 
following types of vacuum integrals: 
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! R (R 2 r[(R-P) 2 ]n 

In dimensional regularization these can all be related to I m n , which in turn reads 



2\-Z-m-n 



r(m + n - f )r(f - m)r(f - n 



(47T 



T(D — m — n)T(m)T{n) 



(3.6) 
(3.7) 
(3.8) 

(3.9) 



2 Various recipes for this can be found in the literature, but we have established all relations from scratch. 
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2 

The remaining P-dependence often appears in the form [(P ■ Q) 2 ]q = ^(gz^j ~~ p»0> where 
we made use of rotational symmetry. 

Within the last part of eq. (|3.4p . a similar expansion can be carried out with respect to 
R in l/(R — P) 2 , etc, but of course not in l/(Q — R) 2 , since both variables are of the same 
magnitude in this case. The challenge then is to deal with the angular integrals over the 
directions of q, r. Letting z = q-r/(gr), q n = a iq, r n = pir, a = ±, p = ±, the most difficult 
structures are the ones containing 



1 



(Q - R) 2 2qr(ap - z) 
The averaging denoted by [...]q,h leads to 



(3.10) 



a=± p=± 

7 E E = 7 E E = . ( 3 - 12 ) 

cr=± p=± (T=± p=± r 

iEE^TT = T~2 • (3-13) 

cr=± p=± 

Antisymmetry kills the structure linear in z, unless there is another angular variable appearing 
in the numerator. Fixing the directions of p, r and integrating over those of q, the latter 
case produces 

where we made use of rotational symmetry and the dimensionally regularized angular inte- 
gration measure. A subsequent integral over the directions of r might contain 

The case without any angular variable in the numerator yields 

1 \ D -3 



1) . (3.16) 



1 - z 2 j D — 4 

Setting finally D = 4— 2e, we obtain the expansions listed in appendix A (eqs. (|A.15|) - (|A.32|) ), 
4. Euclidean momentum-space correlators 



Inserting into eqs. (|3.ip . (|3.2p the asymptotic expansions of the master sum- integrals from 
appendix A as well as the bare gauge coupling from eq. (j2.5j) : expanding in e; and inserting 
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1 = fj,- 2e n 2e " j we can write the results in terms of the renormalized coupling as follows: 
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(4.1) 



(4.2) 



Here and in the following, terms of 0(e) have been omitted where irrelevant. 

The vacuum parts, i.e. the first rows of eqs. ()4.ip . (I4.2p . can be compared with ref. |13j . If 
we rewrite (ji/P) a = (fie/P) a [l — «+...] and re-expand the square brackets, we can reproduce 
the results of ref. [13] . 

As far as the thermal parts go we note that, within the accuracy of our computation, they 
can be expressed as 



AG e (P) 



AG X (P) 
-16c 2 5 4 



+ 



p 2 V 3 
2 



pI 



(4tt) 2 V 3 P 2 



2 203 
+ 18 



g 2 h 



1 + S &o In 

2 



, P 



CeP 2 . 

g 2 N c (22 



(e-3p)(T) + 0(g 



1 + 



(An) 



, 347 

y ln p^ + i^ 



(e+p)(T) 

P 2 ) ' 
(e + p)(T) 



C X P 2 



(4.3) 
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Here we have identified the structures of eqs. (|B.2[) . (|B.3|) in the result, realizing that the 
expressions proportional to (^ — p^)/P 2 couple to the traceless part of the energy- momentum 
tensor, T^ u , satisfying (T^ u ) = (Tqq}(5q^5o u — Si^S^/S), whereas the other terms couple to 
the trace part. In addition we have used renormalization group invariance of (e — 3p)(T) as 
well as the theoretical expectation that the Wilson coefficients should be independent of the 
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"soft scale", T, to provisionally add to the results the logarithmic terms on the second rows; 
the coefficients C, x next to P 2 remain undetermined, because fixing them would necessitate 
a 3-loop computation. The qualitative structure of eq. f|4.3[) agrees with that put forward 
in ref. |10j . (Note that the term proportional to e + p does not appear in classic vacuum 
studies like ref. |16j . because it breaks Lorentz symmetry. In fact, expressed in another way, 
(e +p)(T) = Ts(T), where s denotes the entropy density.) 



Short distances 



As a first concrete application, we consider equal-time correlators in configuration space, 
measured recently in ref. [11]. This means that we need to inverse Fourier transform the 
correlators in eq. f|2.8[) in order to get back to the correlators of eq. (|2.7p . Since we are 
applying the inverse transform not to the full result but to an asymptotic expansion valid in 
the regime P 2 3> (irT) 2 , the inverse transform can (and must) be taken at zero temperature 
(i.e. omitting terms like n B (p)). The master formula for an inverse Fourier transform in 
dimensional regularization reads 
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(5.1) 



We recall from eq. (|2.1Up that r denotes a spatial separation in radial coordinates, r = |x|, 
whereas the temporal separation is chosen to vanish as in ref. [TT], r = 0. 

Employing eq. (|5.ip . we obtain for the structures appearing in the vacuum parts (the first 
rows of eqs. (|4.ip . (|4.2p ) the expansions 
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Note that these are proportional to e, which is why we did not need to show terms of 0(1) 
in eqs. (|4.ip . (14. 2p . The structures multiplying e + p in eqs. (14. 3p . (14. 4|) yield 
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(To arrive at these it is helpful to consider the covariant integral f p e lP ' x P^P v /[P 2 ] a first.) 
Finally, the structures multiplying e — 3p in eqs. (|4.3p . (|4.4p contain kind of an ambiguity 
because, taken literally, the terms independent of P yield 5^ D \x), whereas according to 
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eq. (|5.ip we get nothing. This ambiguity is of little significance, however, since the contact 
terms are of no interest at r 7^ 0. In contrast, the terms with ln(/2 2 / P 2 ) as an integrand yield 
a physical behaviour in dimensional regularization: 
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Inserting eqs. (|5.2p - (|5.6p into the inverse Fourier transforms of eqs. ()4.ip . (|4.2p . we get 
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(5.9) 
(5.10) 
(5.11) 
(5.12) 
(5.13) 
(5.14) 



The results of eqs. (|5.7p . (|5.8p can be compared with those given in refs. [TT] . Inserting 
eg from eq. (|2.3j) . our leading order results for the vacuum and trace anomaly parts of Gq 
agree with ref. [9j. In contrast, we find a coefficient of the term proportional to e +p to be 
larger by a factor 2 than in ref. [9j|§ We remark that this coefficient appears both at 0(g 4: ) 
and 0(g 6 ), and we get contributions at both orders which combine to produce the correct 
renormalization group invariant structure. In addition, the leading order result can even be 
worked out exactly: omitting terms that vanish in dimensional regularization, we get 
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^However the correct coefficient (in momentum space) can be found in the more recent ref. [17] . 
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Figure 2: Numerical estimates of the Wilson coefficients 70.j, le-,e+pi 7x;i' 7x; e + p> CI - e 1 s - (EU), 
(|5.10p . (|5.12|) . (|5.13|) . respectively. Thick lines correspond to choosing \x such that the next-to- leading 
order correction vanishes (we refer to this value as fiopt)] thin lines correspond to /I = 0.5/2 O pt or 
p = 2.0/2opt (the dependence on fl is non-monotonic). The gauge coupling is solved from the 2-loop 
rcnormalization group equation, and A^jg = linip_ ! . 00 /2[&o<? 2 ] bl ^ 2b ° exp[— 55^7] • Particularly in the 
%-channel reasonable apparent convergence is observed only at extremely short distances. 



where A (re) is the scalar propagator in coordinate space, 



A(x) = 



n iPx 



pi 
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Airr 



Re 



coth(7rT(r + it)) + 0(e) . 



(5.16) 



Carrying out the derivatives in radial coordinates [d^Ad^A = d^.Ad^A + 2d T d r Ad T d r A + 
(2/r 2 )d r Ad r A+d 2 Ad 2 A] and setting r = in the end, we obtain an elementary if complicated 
expression, whose expansion in a small rT yields 
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4d 
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. (We 



Identifying e+p from eqs. (jB.2|) . ()B.5p we reproduce the two leading terms of eq. ([5T 
can also read from here that the dimensionless expansion parameter is ~ (rirT) 2 .) 

Concerning G x , it is stated in ref. |llj that it has the same Operator Production Expansion 
as Gg, whereas we find a different sign for the term proportional to the trace anomaly, 
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e — 3p. Our finding is consistent with ref. [16] (cf. eq. (3.14) there): if the vacuum terms 
are normalized to be equal in the two cases, as has been done in eqs. (|5.7p . (|5.8p . then the 
trace anomaly terms come with opposite signs, and cancel in the sum. We note, furthermore, 
that the numerical results of ref. [11] show a very different short-distance behaviour of the 
thermal parts of the two correlators (cf. fig. 5), which could also be a reflection of the fact 
that —G x gets a positive contribution from the trace anomaly at short distances whereas Gq 
gets a negative contribution. (It must be stressed, though, that if the Wilson coefficients of 
the structures proportional to e + p and e — 3p were the same, as is the case at leading order, 
then the positive contribution from e — 3p could not overcome the negative contribution from 
e +p in — G x , given that — (e + p) + (e — 3p) = — 4p < 0.) 

We conclude by noting that the Wilson coefficients in eqs. (|5.9p . (|5.10p . (|5.12p . (|5. 13j) go 
beyond the accuracy of the analysis in refs. [91 [TT]. With next-to-leading order corrections 
available, we should be in a position to estimate the coefficients numerically; this has been 
attempted in fig. El Any sort of apparent convergence is only observed at very short distances, 
particularly in the ^-channel; this unfortunate fact may not be totally unexpected [18j . In any 
case, it can be seen that the various Wilson coefficients could be numerically quite different 
even though they agree at leading order. In principle it would be nice to also determine the 
Wilson coefficients related to the trace anomaly terms (eqs. (|5.1ip . (I5.14p ) at next-to- leading 
order but this would require a 3- or 4-loop computation. (Since the operator yielding e — 3p is 
Lorentz invariant and the Wilson coefficients should be T-independent, it might be possible to 
extract these coefficients from purely vacuum computations, however we have not unearthed 
literature where this would have been achieved.) 

6. Large frequencies 

We next turn to an "opposite" limit, that of large frequencies but vanishing spatial momenta: 
P = (p n ,0). The results of eqs. (|4.ip - (|4.4p remain valid as a starting point. Furthermore, 
it has been argued in ref. [10J that the asymptotic expansions can be analytically continued 
to Minkowski signature, p n — > —i[uj + i0 + ], extracting thereby the spectral functions, p(u) = 
ImG(— + «0 + ]), even though this requires crossing the light-cone, P 2 = 0, on which 
the asymptotic expansions are certainly not valid. In any case, within perturbation theory, 
the procedure should surely be justified, since we could imagine carrying out the analytic 
continuation before the asymptotic expansion in each individual master sum-integral, and 
expanding only subsequently in a large uj/ttT. 

Let us stress that we assume the analytic continuation to be carried out in the presence of 
an ultraviolet regulator for spatial momenta. The regulator is removed (i.e. the limit e — > 
is taken) only after the analytic continuation. 
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With these qualifications the structures in eqs. (|4.ip - (|4.4p yield (at p = and oj > 0) 
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Thereby we get the spectral functions 
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(6.7) 
(6.8) 
(6.9) 
(6.10) 
(6.11) 
(6.12) 



After inserting eg from eq. (I2.3P the leading-order term of agrees with ref. [,19j, but 
we can now add to that result the first correction. Similarly, the leading-order term of 
7 X; i agrees with a result given in ref. [11], but we can add a correction. The coefficients 
70 ;e + p,70 ;e _3p agree with ref. [10], if the coefficient C introduced there is set to C = 1. The 
results for jq. e + p , jq. e _ 3p have more recently been reproduced in ref. [T7] . We are not aware 
of analogous results in the literature for 7 X; e + p , 7 x;e _ 3 P . 

The coefficients 7x;i> f° r which next-to-leading order values are given in eqs. (|6.7p . 
(|6.10p . are estimated numerically in fig. [3j Like in sec. [5j reasonable apparent convergence is 
observed only extremely deep into the ultraviolet regime. 
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Figure 3: Numerical estimates of the Wilson coefficients jg. 1 , 7 X; i, cf. eqs. (|6.7[) , (|6.fOI) . respectively. 
Unspecified conventions are as in fig. [2] Reasonable apparent convergence is observed only at extremely 
large frequencies, u > 50 (79. t ) or w > 80 A^g (j x - 1 ). 



We note that 7e ;e - 3p> 7x;e - 3p (eqs. (|6.9|) . (|6.12|) ) contain second and higher powers of g 2 , 
so that they should gradually decrease with u as dictated by asymptotic freedom (in analogy 
with fig. [3D. On the other hand, if similar results were available in the "shear channel", i.e. 
for 2-point correlators of the traceless part of the energy-momentum tensor, then the terms 
proportional to e — 3p would allow to fix an unknown constant, denoted by D in ref. |10j . 
which is argued to be "saturated" at that level [10]. In that case there is no g 4 multiplying 
it, so the result is a constant. This implies that such terms require a special treatment in the 
context of sum rules, cf. sec. The importance of determining such terms has recently been 
pointed out also in the context of a lattice investigation [20J. Unfortunately, determining 
the 2-point correlators in the shear channel is technically more demanding than the present 
analysis because of the more cumbersome Lorentz structures appearing in the numerator. 



7. On sum rules 



Starting from the spectral representation of a Euclidean correlator, 

dej p(oj) 



G(p n ) 



7T UJ — tp T: 



(7.1) 
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which is generally valid in the presence of an ultraviolet regulator for spatial momentap and 
setting p n — > 0, we obtain 



Adding the dependence on spatial coordinates (which were suppressed above), we thus obtain 
a "sum rule" , relating an integral over the spectral function to a Euclidean "susceptibility" . 

Now, in practice, various complications can arise in the application of sum rules. One of 
them is of an infrared type, and particularly relevant for the "bulk channel" (our correlator 
Ge): indeed there must be a term oc uj 8{oS) in pg because, if we set p — > from the outset, the 
operator couples to a conserved charge J* Too- Such a term does yield a contribution to both 
sides of the sum rule, however not to the transport coefficient of interest (lim a; _ >0 + pg(u})/u), 
so it would be wise to subtract it. In ref. [T7], this was achieved by defining a modified spectral 
function Ap* which has no oo 5(oj) but nevertheless yields the same transport coefficient. If 
the problematic term is subtracted from the right-hand side, it must also be subtracted from 
the left-hand side, and this yields a modified sum rule (cf. appendix C of ref. |21j). However 
this extra contribution is of 0(g 8 ) and thus not directly visible in our 0(g e ) result. 

Another possible complication, of an ultraviolet type, arises if we wish to remove the ultra- 
violet regularization before applying the sum rule. For instance, the vacuum part then grows 
aspocw 4 sign(w), and both sides of eq. (I7.2p are ill-defined. The problem becomes less severe 
if the vacuum parts are subtracted from both sides of eq. (17. 2p . imposing a relation between 
AG and Ap instead. However, even such a subtraction may not be sufficient to remove a 
constant part, as is argued to be the case in the shear channel [211 HDJ; then eq. (17. ip needs 
to contain a supplementary "contact term" on the right-hand side. 

We now return to the case at hand, where only the former problem should be relevant |1Q|, 
ITT] . The left-hand side of the sum rule reads G$(0), G x (0) in our notation. Inspecting the 
master sum-integrals in eqs. ([ATB - (lATT3ji . most are seen to vanish in the limit P — > 0; in 
fact only X a and X\ give a contribution, and the latter can be reduced to the former through 
eq. (IA.14j) . Inserting these values into eqs. (13. lj) . (13.21) we obtain 



where in the final step we made use of eqs. (|2.3p . (|A.23p . Equation (|7.3p agrees with eq. (|B.4p 
and therefore confirms the sum rule of ref. [22], obtained also directly in lattice regulariza- 
tion [23], up to 0(g e ). Unfortunately, as mentioned, the infrared ambiguities discussed in 
refs. |21[ [17] are of 0(g 8 ) and thus beyond the resolution of our computation. 

4 Equation (|7.1|) can be derived for instance by Fourier transforming the known non-perturbative relation 
in eq. Q1.1D . which is unproblematic if G(r) does not diverge at small distances. 




(7.2) 




(7.3) 



G x (0) = 



(7.4) 



14 



8. Summary and outlook 



Following ref. [llj and other recent works, we have analyzed the ultraviolet asymptotics of 2- 
point correlation functions of the trace of the energy-momentum tensor and of the topological 
charge density in finite-temperature Yang-Mills theory. The tool has been dimensionally reg- 
ularized perturbation theory, pursued up to 2-loop order. We have considered both Euclidean 
(operators separated by a purely spatial separation; sec. [5]) and Minkowskian correlators (op- 
erators separated by a Euclidean timelike separation, subsequently Fourier transformed and 
analytically continued; sec. [6|). 

Through our analysis, we have confirmed a number of expressions in recent literature, 
but simultaneously also refined the determination of the corresponding Wilson coefficients, 
particularly in the Euclidean domain (eqs. (|5.9p . (|5,10p . (|5.12p . (|5.13p ). In this case, we have 
also identified some inaccuracies in the literature; on the qualitative level, our expressions 
appear to be in a somewhat better accordance with lattice data than the earlier ones, because 
they show a qualitative difference between the two channels, with the terms proportional to 
the thermal part of the trace anomaly, (e — 3p)(T), coming with opposite signs. Unfortunately 
the perturbative series for the Wilson coefficients show slow convergence and it appears 
uncertain whether, even with the correct values inserted, they can lead to a quantitative 
agreement with lattice data. (The problem is worse in the case of the topological charge 
density correlator, cf. fig. [2) 

As an outlook, we believe that our results can be refined in a number of ways, some of them 
with reasonable prospects for phenomenological success. First of all, it would be interesting 
to determine the full r-dependence of the correlators up to distances r ~ 1/(ttT) rather than 
the asymptotic expansions at r <C 1/(ttT) as in the present study. This would increase the 
range of applicability of the results and might simultaneously improve on the convergence of 
the weak-coupling expansion, both aspects facilitating a comparison with lattice data a la 
ref. [11]. However, we would suggest carrying out the computation after averaging over the 
time separation, r, rather than setting it to zero, because this makes the correlators more 
analogous to the ones encountered in the context of heavy quarkonium physics [7] , a problem 
of actual phenomenological significance. In any case, eqs. (13. ip . (I3.2|) can serve as starting 
points for such investigations; furthermore, the asymptotic results for the time-averaged 
correlators can still be obtained from eqs. (I4.ip - (l4.4p . simply by replacing the inverse Fourier 
transforms in eq. (|5.ip by three-dimensional ones. 

Second, it would be interesting to determine the full w-dependence of the spectral functions 
down to frequencies u ~ irT, rather than carrying out an asymptotic expansion at to ^> irT as 
in the present study because, as mentioned around eq. (jl.ip . such information may be helpful 
in connection with lattice extractions of the corresponding transport coefficients. The full 
frequency dependence is believed to be quite non-trivial and display e.g. a large cancellation 
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of terms of 0(g 4 T 4 ) [21 [TO]. Once more, eqs. (|3.ip . (|3.2p can serve as starting points for 
this exercise; moreover the asymptotic expansions in eqs. (|6,5p . ()6.6|) could provide for useful 
crosschecks on the final results. 

Third, for theoretical completeness, it would be interesting to determine the Wilson coef- 
ficients more accurately, particularly the missing 0(g e ) terms in eqs. (|5.1ip . (|5.14p . because 
this would allow us to estimate the renormalization scale relevant for the gauge coupling 
multiplying the trace anomaly. (Nevertheless, because of the reservations mentioned in the 
present study in the case of the other Wilson coefficients, it is not clear to us whether this 
could lead to a practically successful comparison with the lattice data of ref. [H].) Similarly, 
of course, it would in principle be nice to know all next-to-leading order corrections in the 
Minkowskian domain, relevant for eqs. (16. 8p . (|6.9p . (I6.11j) . (|6.12p . All of these challenges 
necessitate at least 3-loop computations. 

Fourth, it would be interesting to determine the dependence of the various Wilson coef- 
ficients on the number of massless quark flavours, Nf, not least because this offers a fur- 
ther crosscheck on the applicability of the Operator Production Expansion framework in the 
finite-temperature context. Partial results (including only "sea" quarks), and the problems 
encountered in this computation, are discussed in appendix C. 

Finally, it would be interesting to repeat the study in the "shear" channel, i.e. for the 
traceless part of the energy-momentum tensor. The general techniques of sec. EJ as well as 
some of the actual sum-integrals and asymptotic expansions encountered in the present study, 
might play a role in that context as well. 

We hope to return to some of these topics in future work. 
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Appendix A. Bosonic master sum-integrals 



The bosonic "master" sum-integrals appearing in our computation are defined as 
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In fact there is some redundancy here, because the following relation can easily be established 
through changes of integration variables: 

Xi=X a + X c -X f + X r . (A.14) 

Expanding in a small (vrT) 2 / P 2 as explained in sec. El the sum-integrals in eqs. (lA.ip . 
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(|A.2p can be expressed as 
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we can write the expansions of the next-to-leading order sum-integrals as 
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2-e e(2 - e)P 2 \3 - 2e 

Both structures proportional to P 2 , namely 1S3 appearing in almost every master sum- integral 
as well as 1S4 appearing in I c and Id, disappear from the final result for any e. A useful relation 
allowing to change the basis in the terms proportional to 56 is 
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Appendix B. Basic thermodynamic functions 

We recall here perturbative expressions for a number of thermodynamic potentials needed in 
our analysis. The results are needed up to order 0(g e ) in some cases and can be extracted 
from the explicit 0(g 6 ) results in ref. [21]; however, they could also be deduced through 
renormalization scale independence arguments already from the classic 0(g 2 ) results for the 
thermodynamic pressure given in refs. \25\ 126] . 

As explained in the main body of the text, an essential role in the present study is played by 
the energy-momentum tensor of the thermalized system, which in the plasma rest frame takes 
the form diag(e, — p, — p, — p). We separate this into a traceless part, whose 00-component 
reads 

e-i( e _3p) = |(e+p) ! (B.l) 

as well as a trace part, e — 3p. At next-to-leading order, the combination appearing in the 
traceless part can be written as 
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with d,A = N 2 — 1, whereas the leading-order expression for the trace anomaly reads 

(e - 3p)(T) = 2d A g%N c [ " s(g) ^ , (B.3) 

J Q ,r Q r 



with 60 as defined in eq. (|2.6p . Finally, the temperature dependence of the trace anomaly 
reads 



q,r 
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In the expressions above, n B denotes the Bose distribution, n B (q) = \j{eP q — 1), and the 
integrals can be carried out, 

but it is convenient to be able to recognize the integral representations as well. 



Appendix C. On fermionic effects 



We add here the contribution of iVf massless "sea" quarks to the previous results for the trace 
anomaly and topological charge density correlators. By sea quarks we refer to effects origi- 
nating from the fermionic contribution to the gluon self-energy; "valence" quarks would refer 

to the quark part of the energy- momentum tensor operator. It is often said that the quark 

_ _ ^ y 

contribution to the trace anomaly operator, which is proportional to tp^^D^tp or ip^^ D nip, 
vanishes in the chiral limit thanks to the equations of motion; in topological susceptibil- 
ity there should be no valence quark contribution to start with. Below we find practically 
identical expressions for the two channels, however also structures in the Lorentz-symmetry 
violating part of the result which do not fit the form expected from the Operator Product 
Expansion. This could indicate that there is, after all, some mixing taking place and a non- 
zero valence quark contribution to be added at non-zero T; however, given the fair amount 
of work involved, we have not carried out a systematic investigation to clear up the issue. 

With these reservations, the fermionic contributions to the correlators in eq. (|2.8p read 
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Here the following new master sum-integrals have been introduced: 
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The sum-integral denoted by ${q} goes over fermionic Matsubara momenta. As indicated 
this set is somewhat overcomplete (at least in the absence of a chemical potential, as we have 
assumed to be the case throughout). 

We note that, like in the bosonic case, most of the master sum-integrals vanish in the limit 
P — > 0, relevant for the sum rule (cf. eqs. (17. 3D . (17. 4D ). Non-zero contributions arise only from 
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X a and Z a . In 5G X (0) these cancel; for SGg(0) we obtain (inserting also eg from eq. (|2.3p 



5G e (0) = 8d A g%N i (21 a -l a ) 
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where we made use of f {Q) l/Q 2 = -T 2 /24 and f R l/R 2 = T 2 /12. This result agrees with 
the leading fermionic contribution to T 5 d[(e — 3p)/T 4 ]/dT, and therefore conforms with the 
sum rule of ref. [22] at 0(g e ). 

The large-momentum expansions of the master sum-integrals can be worked out like in the 
bosonic case (cf. sec. E}. Supplementing the structures in eqs. (lA.17p - (lA.22j) with 
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where n F denotes the Fermi distribution, n F (q) = l/(e l3q + 1), we obtain the expansions 
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All structures proportional to P 2 , namely £3,53 and £4, cancel in the final results for any e. 

Inserting these expansions into eqs. (jC.ip . (|C.2p and re-expressing the bare gauge cou- 
pling in terms of the renormalized one according to eq. (|2.5p . we obtain results analogous 
to eqs. (|4.ip . (|4.2p . Alas, many different thermal distributions appear and it is not easy to 
express the result in a concise way. We choose rather to insert explicit values, 
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Thereby the fermionic contributions to eqs. (|4.1|) . (|4.2p become 
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Inspecting the results, the terms on the last rows of eqs. ()C,36p . (|C.37p are easy to under- 
stand: the leading-order quark contribution to the trace anomaly is 



e-3p\ 5d A g 4 b N { 



rp4 



288 



(C.38) 



and these terms amount to ^f25(e — 3p)/dAg 2 bo, in perfect accordance with the bosonic results 
of eqs. (|4.3p . (|4.4p . In contrast, the terms on the second rows of eqs. (|C36[) . ()C.37p fit no 
simple pattern. In principle we might expect a fermionic contribution to the Wilson coefficient 
multiplying the leading bosonic e +p, viz. 8ciA7r 2 T 4 /90, as well as direct fermionic effects to 

'e + p\ _ 7N c N { tt 2 5d A g 2 N { 

ua~ ' (C - 39) 

the former multiplied by the next-to-leading order Wilson coefficient and the latter by the 
leading order one. However, there is no way to reproduce the leading fermionic contribution 
to e + p (the first term in eq. (|C.39P ) from the effects in eqs. (|C.36p . ()C.37p (powers of g 2 
and/or group theory factors do not match) and, conversely, there is no way to understand the 
appearance of the temperature-dependent logarithms in eqs. (|C.36p . ()C.37p in terms of the 
Operator Product Expansion contributions. Something is clearly missing and, as mentioned 
at the beginning, one possibility could be an unaccounted mixing with fermionic operators. 
Let us stress that the problem only appears in the contributions proportional to e + p, which 
vanish at zero temperature due to Lorentz symmetry (cf. the last lines of sec. HJ. 
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